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The relationship between natural orbitals, one-body coherences and two-body correlations is ex¬ 
plored for bosonic many-body systems of definite parity with two occupied single-particle states. 

We show that the strength of local two-body correlations at the parity-symmetry center character¬ 
izes the number state distribution and controls the structure of non-local two-body correlations. A 
recipe for the experimental reconstruction of the natural orbital densities and quantum depletion 
is derived. These insights into the structure of the many-body wave-function are applied to the 
predicted quantum-fluctuations induced decay of dark solitons. 

PACS numbers: 03.75.Lm, 67.85.Be, 67.85.De 


I. INTRODUCTION 

In an ideal Bose-Einstein condensate, all bosons oc¬ 
cupy the same single-particle state 4>o{v), whose density 
can directly be inferred from an absorption imara mea¬ 
surement of the reduced one-body density /9i(r) Q. Yet 
in a non-ideal world, not only do interactions between the 
atoms affect the shape of the condensate wave-function 
(()o(r) but also bring more single-particle orbitals into 
play, even at zero temperature, such that pi(r) equals 
an incoherent superposition of their densities in general. 
Theoretically, the many-body state can be characterized 
by the natural orbitals (NOs) [^, i.e. eigenvectors 

of the reduced one-body density operator pi, and their 
populations, i.e. the corresponding eigenvalues: Given a 
sufficiently large weight, the NO of largest population is 
identified with the condensate wave-function and quan¬ 
tum depletion manifests itself in the population of other 
NOs Q. As a matter of fact, correlation effects can be 
traced back to both the occupation number distribution 
of the NOs and their spatial shape, allowing for a mi¬ 
croscopic understanding of various phenomena. Exam¬ 
ples for such phenomena are the Mott-insulating phase, 
where e.g. an intra-well Tonks-Girardeau transition for a 
filling factor of two can be understood in the NO frame¬ 
work Bi, the quantum-fluctuations induced decay of 
dark solitons @-0, where a NO of particular shape is 
dominantly responsible for the soliton contrast reduction 
in the reduced one-body density, and fragmented conden¬ 
sates [lol - [T^ , which are defined as many-body states with 
two or more macroscopically occupied NOs. While there 
are pro pos als for the detection of fragmentation and its 
degree (l^ . [IB |. the one-body density pi(r) has, to the 
best of our knowledge, not yet been unraveled into the 
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contributions |^i(r)p of the individual NOs by means of 
a measurement protocol. 

In principle, the NOs can be obtained from a tomo¬ 
graphic reconstruction of the reduced one-body density 
matrix pi (r, r') and diagonalization. Such quantum state 
tomography is a well established technique for qubit sys¬ 
tems and quantized light fields (l^ . [T^ . For interacting 
ensembles of ultracold atoms, information about the off- 
diagonal elements pi(r,r') can be extracted experimen¬ 
tally from the contrast of interfering slices coupled out 
a trapped Bose gas [isl - f^ or homodyning in uniform 
systems, where one generates copies of one’s system by 
Bragg pulses and let them interfere [2l|, [U. Further¬ 
more, there are theoretical proposals for the pi(r,r') re¬ 
construction based on Raman pulse sequences [23| , vari¬ 
able time-of-flight (ToF) measurements a, heterodyn¬ 
ing with an auxiliary Bose-Einstein condensate [1^ or 
Jaynes principle of maximum entropy [2^. Yet due to 
its non-local character, it is notoriously difficult to in¬ 
fer pi(r,r') experimentally, in particular for non-uniform 
systems. 

Rather than aiming at a completely general reconstruc¬ 
tion scheme for the NOs, we focus here on bosonic many- 
body systems of definite parity with two occupied single¬ 
particle states. Assuming only two occupied orbitals con¬ 
stitutes the simplest, natural extension for bosons beyond 
the mean-field approximation and is physically justified 
in various situations. For this class of systems, we de¬ 
rive an experimentally accessible reconstruction recipe, in 
which density-fluctuation measurements play a decisive 
role, and also gain insights into generic properties of two- 
body correlations. In particular, we show how the char¬ 
acter of the number state distribution function, the NO 
densities, the structure of two-body correlations and the 
relationship between one-body coherences and non-local 
two-body correlations crucially depend on the strength of 
two-body correlations at the parity-symmetry center. All 
these relations are derived exactly from the structure of 
the many-body wave-function. In order to show the im- 
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portance of our results as well as their validity when fur¬ 
ther NOs are slightly populated, we apply our analytical 
methodology to the analysis of numerical ab-initio data of 
the quantum-fluctuations induced decay of dark solitons 
obtained by the Multi-Configuration Time-Dependent 
Hartree Method for Bosons (MCTDHB) (27l - [^ and ad¬ 
dress also the impact of finite experimental resolution. 

This work is organized as follows: First, the setup is in¬ 
troduced in Sec. m Thereafter, we discuss the NO decom¬ 
position of the reduced two-body density in Sec. ImAl 
which forms the basis for the reconstruction of the odd 
and even parity NO density as well as for the characteri¬ 
zation of the spatial structure of two-body correlations in 
Sections IIIIB1IIII Cl and IIIIDl respectively. Finally, we 
apply our insights to decaying dark solitons in Sec. IIII El 
and conclude with Sec. IlYl 


II. SETUP 

In the following, we consider a system of N bosons 
which are energetically or dynamically restricted to oc¬ 
cupy only two single-particle states of opposite parity, 
n\(j)i) = (—l)*|^i), i = 0,1. Here, tt denotes the single¬ 
particle parity operator which inverts either all coordi¬ 
nates or only x. For simplicity, we suppress the y, z ar¬ 
guments in the position representation and remark that 
our results are valid for one-, two- and three-dimensional 
quantum gases. The N^-body state is assumed to pos¬ 
sess a definite A^-body parity, 7rr|d>) = n|d>), H € 

{ — 1,1} with TTr acting on the r-th atom, as it is the case 
for a non-degenerate ground state of a parity-symmetric 
many-body Hamiltonian e.g. of a bosonic Josephson junc¬ 
tion [13 ■ For this class of systems, the many-body state 
is of the form 


K 

\'i') = '^Ak\no{k),N - no{k)), ( 1 ) 

k=0 

where |no,ni) denotes a number state with Ui bosons in 
\(j)i). In the cases of N even and H = 1 or odd and 
H = — 1, the correct parity is ensured by UQ^k) = 2k. 
Otherwise, nQ{k) = 2k-\-l has to be chosen. In all cases, 
K denotes the largest integer with nQ{K) < N. By trac¬ 
ing out A^ — 1 bosons, one obtains for the reduced one- 
body density operator pi = A|^o)(</>o| + (1 ~ 
so that the NOs are given by {(pi). Here, we have in¬ 
troduced the average fra ction of bosons in the even or¬ 
bital, A = rio/N, where (...) denotes the average w.r.t 
the number state probability distribution Thus, 

the quantum depletion equals min{A,l — A} and the 
reduced one-body density is given by the incoherent su¬ 
perposition 

pi{x) = {x\pi\x) = A|())o(a;)|^ + (1 - A)|()>i(x)|^. (2) 

By measuring pi{x) and the real-valued off-diagonal 
elements pi{x,—x) = {x\pi \ — x) for all a;, one could in 


principle reconstruct the NO densities and A without 
knowledge about the full density matrix pi(x^x'). As a 
consequence of the NO parities, one finds \pQ/i{x)\^ oc 
Pi{x) ± Pi (a;, —x) and A = [1 -|- f dx pi(x, —a;)]/2 =[l-l- 
tr(7r/3i)]/2, which links A to the average single-particle 
parity. The drawback of this scheme, however, lies in the 
fact that it requires precise knowledge about pi(x,—x) 
for all X, which is a challenging quantity to measure. 


III. RESULTS 


A. NO decomposition of the two-body density 


Since two-body correlations will indeed provide us 
an alternative pathway to the NO reconstruction, we 
investigate here how the structure of the many-body 
wave-function © manifests itself in absorption image 
noise correlations. The latter have theoretically been 
proven to give valuable physical insights in particular 
for low-dimensional systems [dll - l^ and are measurable 
both after ToF and in situ nowadays [37l - l48l |. For this 
purpose, we derive the two-body density ^2(2^1,2^2) = 
{'ip^{xi)^p^{x 2 )i>{x 2 )'ip{xi))/[N{N — 1)], where 'ip{x) de¬ 
notes the bosonic field operator: 

P 2 {xi,X 2 ) = 2ifi{a(pii{xi,X2)(poo{xi,X2)) + (3) 

+P\(pOo{xi,X 2 )\'^ + 2'^\(poi{xi,X2)\'^ + S\(pii{xi,X 2 )\'^, 

with pij{xi^X 2 ) abbreviating the normalized sym- 
metrization of the Hartree product (pi{xi)(pj{x 2 ). In the 
following derivation, we will eliminate the off-diagonal 
term with the coefficient a, which is a function of the co¬ 
herences between the respective number states, 

by virtue of the parity symmetry. The second coeffi¬ 
cient is related to the second moment of the number 
state distribution \Ak\^ via /3 = [rig — no]/[N{N — 1)] 
and determines the remaining coefficients 7 = A — /3 and 
(5= l-h/3-2A. 

Assuming a finite central density, pi(0) > 0, we cal¬ 
culate the two-body correlation function 52(2:1,2:2) = 
52(2:1,2:2)/[51(2:1) Pi (2:2)] (4^. at the symmetry cen¬ 

ter 


52(0,0) 


A 

A2 


N 


N-1 



var(no) - no\ 

J 


(4) 


where var(rig) = ~ ^o- measuring the central 

density and its fluctuations, one can in principle deduce 
52(0,0) and, thereby, characterize the number state dis¬ 
tribution in the categories Poissonian, sub- and super- 
Poissonian. Similarly, a measurement of the n-th order 
correlation function gn{xi = 0, ...,Xn = 0) gives insights 
into the n-th moment of the number state distribution 
(cf. e.g. the experiment El for n = 3). As we will show, 
the strength of two-body correlations at the symmetry 
center constitutes a key parameter, which controls both 
the reconstruction of the NO densities and the relation¬ 
ship between local and non-local two-body correlations. 
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Regarding the impact of an unavoidably given finite ex¬ 
perimental resolution, we refer to the discussion at the 
end of Sec. nTTEi 


B. Reconstruction of odd NO 


In order to eliminate the term proportional to a 
in we make use of (^i(O) = 0 and consider the 
non-local two-body correlations g2{x,0) = [P\(j)o{x)\^ + 
7|^i(x)p]/[Api(x)]. After substituting the density of the 
even NO \(i)o{x)\^ via ([ 2 ]) and employing ([U, we obtain 




32 ( 0 , 0 ) - g2{x,Q) 

32 ( 0 , 0 )-! 


Pi{x), 


( 5 ) 


which hold^ for non-trivial two-body correlations at the 
symmetry center, i.e. 32(0,0) ^ 1 . Under this condition, 
we have thus shown that the density of the odd NO is pro¬ 
portional to the total reduced one-body density spatially 
modulated by the strength of non-local two-body correla¬ 
tions between the symmetry center and the position x of 
interest. This relationship constitutes a key result of this 
work since it provides an explicit reconstruction scheme 
for the microscopic quantity |^i (x) p in terms of the mea¬ 
surable densities pi{x) and p2(x,0). In particular, this 
simple reconstruction recipe does not require to measure 
the off-diagonal elements pi(x, —x). For the ground state 
of a parity symmetric Hamiltonian with short-range in¬ 
teractions, the many-body parity turns out to be even 
(5^ and, for not too strong interactions, most of the 
bosons occupy the NO of even parity. Thus, m gives 
experimental access to the density of the orbital whose 
population is responsible for quantum depletion. Besides, 
the positive semi-definiteness of |^i(x)p and pi{x) im¬ 
plies that (anti-)bunching at the symmetry center yields 
32(0,0) as an upper (lower) bound for the non-local cor¬ 
relation 32(x, 0). 


C. Reconstruction of even NO 

Due to the presence of the unknown A in ( 2 ]), den¬ 
sity and density-fluctuation measurements can only be 
employed to relate the NO density difference |(/)o(x)p — 
1^1 (x)p at different points in space, x = Xi/2, but not to 
extract | 0 o(a^)P itself. Nevertheless, using ® the pos¬ 
sible candidates for |0o(x)p can be restricted to a one- 
parametric famil} 13 , 

\ Mx )\^ _ f 32(0,0)-32 (a^, 0 ) \ 

pi(x) U ^ 32(0,0)-! J ’ 


^ In the limit 92 ( 0 , 0 ) —>■ 1, the right hand side of l[5j tends to 
such that no information about the NO is gained. This 
fact is also reflected by 92 ( 0 , 0 ) = 1 implying 92 ( 0 :, 0 ) = 1 (cf. 
Sec. ImPll . 

^ In the 92 ( 0 , 0) —> 1 limit, this equation becomes equivalent to l[2j 
such that no information is gained. 


where A £ ( 0 , 1 ). Thus, a theoretical estimate for A 
by means of e.g. number-conserving Bogoliubov theory 
| 53 h 571 | would allow for uniquely determining the density 
of the even NO. In order to obtain a measurement proto¬ 
col for A as an alternative, we ins pect the first- order co¬ 
herences 3i(x, —x) = 3i(x, —x)/yjpi{x)pi{—x) (4^. [^. 
Substituting |(/)o/i(x)p by (0 and (O, respectively, in 
3i(x, —x), we arrive at 


3i(x,-x) 


1 0^1 A^ff2(0,0)-32(x,0) 

^ ^ 32 ( 0 , 0 )-! 


( 7 ) 


Thus, the additional knowledge of 3i(x*,—x*) for some 
convenient position x* with 32(x*,0) 7^ 32(0,0) is suffi¬ 
cient for determining A. In the case that we have only 
experimental access to the modulus of 31 (x, —x) but not 
to its sign, we may extract A = A±(x) for both signs, 
i.e. ±|3i(x, —x)|, from ([ 7 ]) for all x of finite density with 
32(x, 0 ) 7^ 32(0,0). One easily verifies A+(x) > A_(x) 
and, in many situations, the local sign of 31 (x, —x) can 
then be fixed by requiring A not to depend on x. Know¬ 
ing A and an estimate for A, we may also infer var(no) 
from ([ 3 ]). Finally, Eq. ([ 7 ]) gives the conceptual insight 
that the average fraction of bosons in the even NO me¬ 
diates a relationship between the first order coherences 
31 (x, —x) and two-body correlations 32(x,0). 


D. Spatial structure of two-body correlations 

While the two-body correlation function features a 
particle exchange and a two-body parity symmetry, it 
does not remain invariant under a parity operation act¬ 
ing on one atom only. By inspecting p2(xi,X2) -I- 
32(xi,—X2) — 23i(xi)3i(x 2), i.e. essentially the sum of 
density-density correlations at (xi,±X2), the parities of 
the NOs can be employed to eliminate the off-diagonal 
term oc a such that a relationship between 32(xi,X2) 
and 32(3^1,—X2) can be established. Here, we have to 
distinguish two cases: (i) In the absence of two-body 
correlations at the symmetry center, i.e. 32(0,0) = !, 
we obtain the relation 32(xi,X2) -I- 32(xi, — X2) = 2 , 
which has three important consequences. Firstly, the 32 
function is fully determined by its values in the sector 
S = {(xi,X2)|0 < xi < X2} (cf. Fig. [T](a) for an illustra¬ 
tion). Secondly, local bunching (anti-bunching) struc¬ 
tures 32 (xi, X2) > ! (32(3^1, X2) < !) forxi ~ X2 translate 
into non-local anti-bunching (bunching) structures of the 
same magnitude at (xi, — X2). Thirdly, pairs of atoms are 
uncorrelated on the Xi/2 axis, i.e. 32 (x, 0) = 32(0, x) = !. 
(ii) In the presence of two-body correlations at the sym¬ 
metry center, we may employ the reconstruction formula 
([ 5 |) to obtain the following functional equation, which has 
to be fulfilled for every 32 with 32(0,0) 7! ! in order to 
be compatible with the many-body wave function m, 

32(3:1, X2)-b 32(2:1,-X2) = 2(^1-b (8) 
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Figure 1. (Color online) (a) Structure of two-body correlations for the case 32(0,0) = 1. (b) g 2 {xi,X 2 ) at t = 5.5rc for a 
density-engineered initial state oi N — 100 bosons with 7 = 0.04 in a one-dimensional box of length L = 20^ (obtained by 
(ML-)MCTDHB simulations). The dashed lines indicate the positions of two counter-propagating gray solitons and the arrows 
point into their directions of motion. The circles indicate the characteristic correlation pattern observed for a single gray soliton 
in 0. Strength of two-body correlations at the symmetry center: 32(0,0) « 1.003. 


with f{x) = g 2 (x, 0 ) — 1 . This restriction on the func¬ 
tional form of 32 may be used experimentally as a neces¬ 
sary condition for testing the validity of the two-orbital 
approximation. 


E. Applications 

Dark solitons, being well-known for their stability 
within the mean-field approximation (see and ref¬ 
erences therein), suffer from a quantum-fluctuations in¬ 
duced decay due to an incoherent scattering of atoms 
from the soliton orbital into an orbital localized at the 
soliton position, c.f. e.g. [M. Due to the increas¬ 
ing population of this orbital, the depth of the char¬ 
acteristic minimum in the reduced one-body density is 
reduced, i.e. the soliton contrast decreases on average 
over many absorption image measurements. Since this 
decay process can qualitatively be understood within 
a two-orbital approximation, dark solitons constitute a 
straightforward example for testing the validity of the 
above insights in situations when further orbitals par¬ 
ticipate with, however, minor weight. In the following, 
we consider N bosons of mass m in a one-dimensional 
box potential of length L with a contact interaction 
strength g. This system is governed by the Hamiltonian 
H = + in a unit system based 

on the chemical potential po = gN/L, the healing length 
^ = h/^mgo and the correlation time Tc = fi/go, where 
7 = mgL/ih^N) denotes the Lieb-Liniger parameter. 

Both for finding the initial state and for the sub¬ 
sequent propagation in the following two scenarios, 
we employ our recently developed Multi-Layer Multi- 
Configuration Time-Dependent Hartree Method for 
Bosons (ML-MCTDHB) [2^, . which reduces to the 

pioneering MCTDHB method [ 23 | if applied to a sin¬ 
gle species in one spatial dimension. This method is 


based on an expansion of the total many-body wave- 
function with respect to bosonic number states with an 
underlying time-dependent, dynamically optimized sin¬ 
gle particle basis of M states. All conceivable number 
state configurations for the given M single-particle states 
are into account. Being based on a variational princi¬ 
ple, the (ML-)MCTDHB equations provide us with a 
variationally optimized solution to the time-dependent 
many-body Schrodinger equation. By incrementing M, 
we found in both scenarios discussed below that conver¬ 
gence is ensured on the considered time-scale if M = 4 
(cf. 0 ). A NO analysis of the full numerical data reveals 
that two NOs contribute with significant weight while the 
probability to find an atom in one of the other two NOs 
does not exceed 0.023 for the time-scales of Fig.[T](b) and 
Fig.H 

First, we start with the ground state oi N = 100 
atoms in the box with an additional Gaussian barrier 
V{x) = hexp[—a:^/( 2 r(;^)], h = 60 po, w ~ 0 . 07 ^ so that 
we engineer a pronounced density notch at x = 0 . Hav¬ 
ing switched off the barrier, we let the interacting many- 
body system evolve in the box potential. In the course 
of time, the single density minimum splits into a pair of 
counter-propagating gray solitons, which are slowly de¬ 
caying due to quantum-fluctuations [^. We have shown 
that a single gray soliton is accompanied by highly local¬ 
ized two-body correlations in the vicinity of the instan¬ 
taneous soliton position resulting in a bunching of atoms 
in the soliton flank opposite to its direction of move¬ 
ment (cf. circles in Fig. [Ijb)). Yet the 32 function of two 
counter-propagating gray solitons turns out to be more 
than the sum of the local correlation patterns of the in- 
dividnal gray solitons - additional non-local correlations 
occur between the two solitons (Fig. [ijb) and Q). Ob¬ 
serving numerically that |32(0,0) — 1 | <C 1 , we may now 
understand these non-local correlations as a generic prop¬ 
erty of parity-symmetric systems with essentially two oc- 
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cupied NOs: Denoting the position of the soliton mov¬ 
ing to the right/left with our above results imply 

g2{xR - e, XL + e) ~ 2 - 32(2;^ - e,XR - e) such that 
bunching in the back of a single soliton (e > 0) translates 
into antibunching of approximately same magnitude for 
finding an atom each in the back of each soliton. 

Second, in order to realize a situation with significant 
deviations of 32(0, 0 ) from unity, we additionally imprint 
a relative phase of tt between the two half-spaces at t = 0 
[M. Thereby, a black soliton is initialized at a; = 0 . 
As time evolves, the density minimum becomes filled up 
by incoherently scattered atoms inducing strong bunch¬ 
ing correlations at the symmetry center. Fig. [ 2 ] clearly 
shows that the density of the dominant NO of odd par¬ 
ity, |(/)i(a;)p, features the characteristic density notch of 
a black soliton and can be reliably reconstructed by the 
scheme m at times when the soliton contrast in the full 
density pi{x) has been reduced. For longer times, the re¬ 
constructed 1^1 (a;) p deviates slightly more from the full 
numerical results since the third and fourth dominant 
NO have gained more population. The reconstruction 
of \(j)o{x)\'^, i.e. the NO mostly responsible for the soli¬ 
ton decay and strong two-body correlations, turns out to 
be more sensitive to the slight population of these NO: 
Due to the phase-imprinting scheme, we expect most of 
the atoms to reside in an odd NO. Finding numerically 
A_(a;) < 1/2 < A+(a;), we thus take the negative sign 
of gi{x,—x) and estimate A by averaging A_(a;) over 
some interval. As a result, we can fairly well recon¬ 
struct |</o(a:)P according to (|6l) for not too long times 
(Fig. Ela)). Thus, our scheme can be used to experi¬ 
mentally verify the microscopic decay mechanism of a 
black soliton contrast in the reduced one-body density 
via a NO being localized at the position of the soliton, 
given that thermal excitations are sufficie ntly suppressed 
as achievable in nowadays experiments | 59 l |. At longer 
times {t > 3 tc, see Fig. E^b)), i.e. when the assumption 
of only two occupied orbitals becomes less valid, however, 
the reconstructed \(j)Q{x)\'^ deviates much stronger from 
the full numerical results compared to the reconstruction 
of \Mx)\^- 

The reconstruction formulas ® and ([ 5 ]) have been de¬ 
rived under the assumption of perfect experimental reso¬ 
lution. In order to test the robustness of the reconstruc¬ 
tion recipe against finite resolution, we have (i) embed¬ 
ded our numerical simulation results for the considered 
purely one-dimensional model into three-dimensional co¬ 
ordinate space by assuming all atoms to reside in the 
ground state of the transverse harmonic oscillator and (ii) 
convoluted the one- and two-body densities with a Gaus¬ 
sian point-spread-function of variable width a. These 
coarse-grained quantities have then been inserted in the 
reconstruction formulas m and (H]). We have found 
that the NO densities and thereby the decay mecha¬ 
nism of the soliton contrast in the reduced one-body den¬ 
sity can be reconstructed as long as a < ^/ 4 , which 
was to be expected (plots not shown). This require¬ 
ment on the optical resolution is demanding, of course. 




X (units of 0 

Figure 2. One-body density and density of the two most 
dominant NOs at times t = 2 tc (a) and t = 3.5rc (b) for a 
many-body system initially featuring a black soliton centered 
at X = 0. All parameters are as in Fig. [Ub). Dotted line: 
pi(x). Solid lines: Exact NO densities |(^o(®)|^ (black) and 
\4>\{x)'^ (gray, reduced by factor 2). All these curves are ob¬ 
tained from (ML-)MCTDHB ab-initio simulations. Dashed 
lines: Corresponding reconstructions of the NO densities. Es¬ 
timated A « 0.079, 0.088 for (a), (b), respectively. Exact 
values: A « 0.073, 0.098. Two-body correlations at the sym¬ 
metry center: 32 ( 0 , 0 ) « 2.699, 2.013. 

but not out of reach since the healing length could be 
raised above the resolution limit by a Feshbach reso¬ 
nance, for example. For such values of cr, the recon¬ 
structed |(/o(a;)P density happens to be sufficiently ro¬ 
bust against errors in the depletion A, for which we 
assumed a relative uncertainty of 20 %. Finally, we re¬ 
mark that not ^2(2:1,2:2) but density-density correlations 
([^/■|■(xl)1/(xl) - Npi{xi)][tp^X2)'il}{x2) - Npi{x2)]) are 
measured in experiments. Due to the canonical com¬ 
mutation relations for bosonic field operators, however, 
these two quantities essentially differ by an autocorre¬ 
lation peak (X ( 5 (xi — X2), which becomes smooth when 
convoluted with a point-spread-function and whose con¬ 
tribution is suppressed as 1 /A for large particle numbers. 

IV. CONCLUSIONS 

We have shown how physical knowledge about the 
structure of the many-body wave-function can be em¬ 
ployed for deriving generic properties of two-body cor¬ 
relations and an experimentally relevant reconstruction 
scheme for the NO densities. In addition to the discussed 
dark soliton example, our results should be applicable 
to many other systems such as Bose gases in a double¬ 
well potential [33, fragmenting bright solitons (T^ . [T^ 
or s ymm etrically colliding fragments in a harmonic trap 
[la [53 ■ If the central density turns out to be vanish¬ 
ing or too small such that 32(0,0) becomes effectively 
ill-defined, the whole analysis has to be carried out in 
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momentum space via long ToF measurements. We hope 
that our work stimulates the interest in NO reconstruc¬ 
tion schemes such that these microscopic quantities be¬ 
come experimentally accessible for a broad class of sys¬ 
tems. 
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